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Abstract 

We determine non-perturbatively the fixed-point action for fermions 
in the two-dimensional U(l) gauge (Schwinger) model. This is done 
by iterating a block spin transformation in the background of non- 
compact gauge field configurations sampled according to the (perfect) 
Gaussian measure. The resulting action has 123 independent cou- 
plings, is bilinear in the Grassmann fields, gauge-invariant by con- 
sidered the compact gauge transporters and localized within a 7 x 7 
lattice centered around one of the fermions. We then simulate the 
model at various values of (5 and compare with results obtained with 
the Wilson fermion action. We find excellent improvement for the 
studied observables (propagators and masses). 
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1 Introduction and Motivation 



It is expected that lattice actions, which lead to continuum theories at crit- 
ical points, obey universality. However, since computer studies near critical 
points are plagued by various obstacles like finite size effects and critical slow- 
ing down, it makes sense to work with lattice actions that reproduce general 
continuum properties in some sense faster. Technically speaking, these ac- 
tions have smaller corrections in powers 0(a n ) of the lattice spacing constant. 
Optimally, such an action has no such corrections and thus no corrections 
to the leading critical behaviour, i.e. no corrections to scaling. Such actions 
have been called "perfect". 

The traditional lattice actions suggested for gauge theories are simply 
confined to a few lattice field variables (are ultra-local), but have corrections 
0(a 2 ) for bosons and even larger corrections 0{a) for fermions. The im- 
proved actions have to introduce more terms. As long as the contributions 
are exponentially damped with regard to their extension in real space one 
calls the action local. More terms complicate the simulation and eventually 
one has to find a compromise between efficiency and perfectness. Most of the 
improvement programs add to traditional actions a few terms adjusting the 
couplings such as to cancel lattice correction in some 0(a n ). This was done 
within lattice perturbation theory PJE] or relying on non-perturbative meth- 
ods for determining the weight parameters for the various terms in the action 
[3]. There are also non-perturbatively motivated suggestions for tadpole- 
improved actions [3]. 

The other group of improvement programs is inspired by scale transfor- 
mations, in particular block spin transformations (BST). Asymptotically free 
theories have their continuum limit at vanishing gauge coupling. In such a 
situation one may identify a fixed point (FP) of a BST by the solution to the 
classical field equations of the combined action and BST of the spin model 
jS]. For quadratic actions the fixed point may be obtained quasi analytically 
(cf. [6] for the Gaussian model and [Zj for free fermions and gauge fields). 
Another elegant approach is direct blocking from the continuum theory jHllH]- 

The lattice actions obtained in this way have the remarkable property that 
they are classically perfect, in the sense that the solutions of the equations 
of motion are related to their continuum counterparts. Furthermore the FP- 
action is tree-level Symanzik-improved to all orders in the lattice spacing a 
[10J and there is evidence, that cut-off effects are strongly reduced also at 
the 1-loop-level HH Q2J EE! . 
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Our approach belongs to the second mentioned type of programs. For the 
2D U(l) gauge theory with Wilson fermions (the lattice Schwinger model) we 
determine non-perturbatively an optimal fermion action in the background 
of gauge field configurations sampled according to their (optimal) Gaussian 
measure. We choose the massless Schwinger model as our testing ground 
since there we have some experience and the possibility to compare with 
other results. FP-actions for that model were also studied by |H1 Clj with the 
method of small fields and recently by ^3] in a perturbative expansion. 

We are interested in some lattice representation of the continuum action 
for the massless model JE] 

X - J d 2 xF(x) 2 - J ^x^ixh^d^ + tgA^^ix) . (1) 

with one or two flavours of fermions. The lattice action should respect the 
basic symmetries like gauge invariance, translational and rotational invari- 
ance, parity symmetry, charge conjugation and the hermitian invariance and 
should have the correct (naive) classical limit. Also one has to take care for 
the fermion doubling problem. Beyond these requirements the form of the 
lattice action is largely arbitrary. 

Here we use a real space BST with a blocking factor of 2 (Sec. 2). We 
show how one can simplify the BST in the limit (3 — > oo, and with this approx- 
imation we iterate the BST to determine a FP-action. The parameters are 
chosen such that the BST yields the most local action in the non-interacting 
case. Since we use the non-compact formulation, the gauge field part of the 
FP-action may be determined analytically jH]- The fermionic part has to 
be determined numerically. The resulting FP-action in this approximation 
defines a classically perfect action and for large (3 we expect that it is a good 
approximation for the renormalized trajectory. In this limit the fermions de- 
couple from the gauge field part of the BST, and the gauge field acts like a 
background field for the fermionic sector JTj, therefore the fermionic action 
stays quadratic in the fermionic field variables. 

We suggest a parameterization for the fermionic action with terms in a 
7x7 square on the lattice, using compact link variables. With this parame- 
terization and the Wilson fermion action as starting point we determine the 
FP-action by iterating the BST. In each step we generate 50 coarse gauge 
field configurations taken from their distribution defined by the gauge field 
FP-action. In the background of these gauge fields we block the fermions 
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from a 14 x 14 lattice down to the 7x7 lattice and determine the iterated 
fermionic couplings. 

In the limit of small A^x) one can compare our FP-action with the one 
obtained by perturbative methods; for the so-called clover term our results 
agree up to a few percent with those published in [T7j. We then simulate 
our improved action at various values of (3 for both, the 1- and the 2-flavour 
Schwinger model (sec. 3). We find substantial improvement of the rota- 
tional invariance of the mesonic 2-point function, continuum-like dispersion 
relations and excellent results for the bound state masses. 

2 Discussion of the Method 

We denote the lattice action by 

/3S{A)-S(V,*,A) , S{$,*,A) = f M(A)tf , (2) 

where S(A) denotes the gauge field part, M the lattice Dirac operator matrix, 
and (3 = 1/g 2 is the gauge field coupling. 

We block from a so-called fine square lattice with sites x G Zn x Zn 
to a coarse lattice organizing the fine lattice in 2 x 2 blocks. These blocks 
constitute the points x' of the coarse lattice. We enumerate the sites x' by 
pairs of odd numbers such that a site x' corresponds to the block (x', x' + 
1, x' + 2, x' + 1 + 2) on the x-lattice. The Grassmann fields are ^(x), \fr(x) 
(respectively W(x'),^'(x') on the coarse lattice); the non-compact gauge 
fields Afj,(x) G M live on the links (x,x + ft). For the fermions we use anti- 
periodic boundary conditions and for the gauge field periodic ones. 

The BST is defined as 

-/3'S'(A')+S'(9',9',A')+c _ 

= J DfADm e ~P(S(A)+T(AA'))+S(*,y,A)+T(V,*^'^') ^ ( 3 ) 

where c is an irrelevant constant and DfA is the measure for the gauge field 
including a local gauge fixing. 

We fix the gauge within each block. A gauge field configuration on the 
fine lattice A^x) is in the so-called fine gauge if and only if 

W : A l (x') = -A l (x' + 2) = -A 2 (x')=A 2 (x' + i) (4) 
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is valid (cf. Fig. QJ. It can be shown by explicit construction that for each 
gauge field configuration A^(x) there exists a unique gauge field configu- 
ration Ap{x) in the fine gauge which is related to the first one through a 
fine gauge transformation. These are gauge transformations, which leave the 
BST-Kernel T(A, A') (for fixed A' , see below) invariant. 




-Ai(k'+2) 
•«« " 



-A 2 (x') 



A 2 (x'+1) 



x=x ' A x {x) x'+l 




Figure 1: In the fine gauge the plaquette field strength within each block is 
distributed equally among the four link variables. 



To get a well defined path integral in the BST and to apply our method 
we have to get rid of these fine gauge degrees of freedom. For that reason 
we integrate not over all gauge field configurations but only over those which 
are in the fine gauge. 

j D l A ==/ DA - (6) 

the fine gauge 

The kernel of the BST for the fermions was taken from |%1 114j. 
r(#, #') = 

£ E (W) - 1 E *«(*)) (*J(*o - 1 E . (6) 

where one has has to choose 62 = v^2 in order to have a fixed point for this 
BST and 02 = 1/4 for maximum locality in the situation of free fermions. 
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Here ^2 x<Ex , denotes the sum over all fine lattice sites x belonging to the 
coarse lattice site x' . 

For the kernel of the gauge field we define an average over the four 2-link 
connections between corresponding sites in adjacent blocks, 

KW = y £(4.(*)+4.(* + A)) ■ (7) 

XdX 1 

A lattice differential operator of second order may be defined as 

A^Bl(x') = V (25>') - B'^x' - 2/2) - B'^x' + 2/2)) . (8) 
With these conventions we write the BST for the gauge fields as follows: 

T (A A') = \Y1 W) - B'^x')) (a 2 + ^A,)- 1 (A^(x') - B'^x')) . (9) 

x',fi 

For the wave function renormalization factor fl 2 , we choose the value 2 (as 
in the compact case of [8j) in order to make the action on the coarse lattice 
gauge invariant. This value is different from the one which would follow 
due to dimensional considerations and, in fact, for (3 2 — 2 there exists only a 
quasi-FP for the BST of the gauge field. For this value of (3 2 for ultra-locality 
one has to choose a 2 = 1/2 and 72 = —1/32. 

With these values the BST respects the basic symmetries of the Schwinger 
model with Wilson fermions. The resulting action on the coarse lattice is 
gauge invariant, hermitian invariant, invariant under the charge conjugation 
and respects the lattice symmetry. It does violate chiral symmetry. 

For given A' there exists a unique minimizing configuration of S(A) + 
T(A, A') which we denote by A m \ n {A ! ). Since we use a non-compact gauge 
field with action and BST quadratic in the fields, A min (A') can be computed 
straightforwardly by solving a set of linear equations. For (3 — > 00 the saddle 
point A min (A') dominates the path integral of the gauge field giving 

e -0'S'(A')+S'^'^',A')+c = e -P(s(A min )+T(A min ,A'))+ci x 

x(^J D'VD^ exp ^, A min ) + T(%, ^, *')] + O 

where c\ is some constant which may be absorbed into c. For Grassmann 
variables the "gaussian" integral results in the exponential of an element of 
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the Grassmann algebra in and This algebra element is a sum of bilinear 
terms in the fermionic variables plus a constant. It may therefore be identi- 
fied with the blocked fermionic action by determining the coefficients of the 
corresponding fermion matrix. Formally this is equivalent to the saddle point 
minimization for the bosonic action. Due to this analogy to bosons we will 
denote this process of integration over the Grassmann fields and subsequent 
identifications of the coefficient of the fermionic action by min.r^ 

Having solved the two path integrals (the bosonic and the formal fermionic 
one) we separate the resulting action of (i.e. the logarithm of ()10j)) into two 
parts, one belonging to the scalar subspace of the Grassmann-algebra, the 
other one quadratic in the fermionic fields. 

P'S'(A') = (3(S(A miQ )+T(A min ,A')) + c 2 (A min ) + 

S'(*',*',A') = min (5(f,*,4 m ) + T(*,f,f',f'))+0^j (11) 

The constant C2(A min ) is just the logarithm of the fermionic determinant 
resulting from the Grassmann integral. Compared to the leading term pro- 
portional to (3 we may neglect this contribution and find that the fermions 
decouple completely from the gauge field BST. The defining equations for 
the FP-action finally have the form 

P'S*(A') = P min^ in the fine gauge (S*(A) + T(A, A')) , (12) 
A') = min { ^ } , A min ) + T(V, , *')) ■ ( 13 ) 

These equations replace now the BST at (3 — > oo. Since we want to use 
the action also at moderate (3 values, we have to calculate S' for strongly 
fluctuating configurations A', too. This however is naturally possible, since 
we never demanded that A' ought to be small. As was already mentioned in 
|Hj this method has nothing to do with perturbation theory. It can be shown 
that the form of the action is independent of the number of fermion species. 

Here we should mention that both, the FP-action for the gauge field and 
the FP-action for the fermions do have scale invariant solutions. For the 
gauge-field this was proven in j^j; for a non-compact gauge theory with pe- 
riodic boundary conditions this is not very interesting, since there are no 
instanton solutions. For the fermionic FP-action the theorem takes the fol- 
lowing form. Suppose that and are solutions of the classical equations 
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of motion dS(W, A')/dW = and dS(W, A')/dW = 0, i.e. are part of 
the null-space of the fermionic matrix, then also the minimizing configura- 
tions of (JT3jl ^f, \& on the fine lattice are solutions of the equations of motion 
corresponding to S(^f, ^, A min (A')). 

In the following two sections we discuss how to solve (fT2*j) by analytical 
methods and (fT3|) by numerical methods. 

2.1 The FP-action for the gauge field 

For details we refer to f3J IE] • For our choice of the parameters of the BST 
the ultra-local standard (non-compact) plaquette action is a fixed point, up 
to the wave function renormalization. In d = 2 this action is 

Sp(A) = ^$>(*) 2 

X 

= l - ^(A 2 (x + l)-A 2 (x)-A 1 (x + 2) + A 1 (x)) 2 . (14) 

x 

As it was shown in [E] on an infinite lattice the equation 

jS P (A') = mm(S P (A)+T(A,A')) (15) 

4 {A} 

is valid. It is remarkable, that (fT5|) is also exactly fulfilled on a finite lattice 
as long as Sp fits on the coarse lattice. For this reason we forget about the 
finite size effects even in the case of the FP-action for the fermions (cf. also 

HE!). 

Although Sp is strictly speaking not a fixed point under the BST (due 
to the necessary but trivial rescaling of (3) it is still a perfect action. In 
the Schwinger model it is classically perfect and for the free gauge field it 
is quantum perfect, since it describes the same physics as on an infinite fine 
lattice. To determine the FP-action for the fermions it is necessary to iterate 
the BST and therefore it is necessary to renormalize (3' — > 4/3 after each step 
to avoid the unwanted convergence of f3 — *■ 0. 

2.2 The FP-action for the fermions 

The fixed point action for the fermions can be calculated in a perturbation 
expansion as was shown in [El E] or more explicitly in But in ^Hl 
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the authors demonstrate for the case of the Gross-Neveu model at large N 
that the non-perturbative approach of |Hj is really necessary to calculate the 
perfect action. 

Eq. ()13|) defines the action on the coarse lattice for any configuration 
A', and - as we discuss later - one can advise numerical procedures which 
calculate the action to high precision. However, since we want to iterate the 
BST and since we want to use this action in numerical simulations we have 
to parameterize S'(j&', ty', A') with a finite number of coupling constants. 

3 

S F (% % A) = VM F (A)y = J2J2 PiW *(x)^U(x, fMx + 5f) . (16) 

i=0 xj 

Here Mp(A) is the parameterized fermion matrix. By / we denote a closed 
loop through x or a path from the lattice site x to x + 5f (Sf denotes the 
distance vector on the lattice corresponding to the path /). U(x,f) is the 
parallel transporter (i.e. the element of the compact group) along this path. 
The o"j-matrices denote the Pauli matrices for i = 1,2,3 and the 2x2 unit 
matrix for i — 0. 

We have to ensure that the BST stays on the critical surface and eventu- 
ally converges to the nontrivial fixed point. Thus we impose the condition 

£>(/) = 0, (17) 
/ 

which guarantees, that Sp reproduces the action of the massless continuum 
Schwinger model in the naive continuum limit. The normalization of p is 
fixed by demanding 

S>(/)0V)i = 1 ■ as) 

/ 

Here (Sf)x denotes the component of the path vector in the 1-direction. We 
want to emphasize that ()16|) is not the most general parameterization. E.g., 
since the gauge fields enter as parallel transporter in their compact version a 
clover-term ^/(x)F fJ-u (x)^/(x), which contains the non-compact field strength, 
can be represented only in an approximate way for values —it < F(x) < it. 
For sufficiently large values of f3 this should be no problem. 

We require the invariance of the action under certain symmetries; thus we 
can impose some conditions for the coupling constants Pi(f). For example 
hermitian invariance and invariance under charge conjugation implies that 
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Po,i,2(/) has to be real and pz(f) has to be purely imaginary. These and fur- 
ther symmetries due to the lattice geometry drastically reduce the number 
of independent coupling constants. We have considered terms that connect 
the central site x with any other site x + 5f in a 7 x 7 lattice. Concerning the 
length of the connecting paths / we first considered paths which may exceed 
the shortest connection between x and x + 5f by up to 8 links and some 
extra paths containing higher powers of plaquette variables. However, in the 
iteration procedure it turned that one may omit many of these. Altogether, 
respecting the mentioned symmetries, we finally considered 33 different geo- 
metric shapes corresponding to 123 independent coupling constants [20J. 

For the determination of the fixed point action we proceed as follows. The 
starting point for the fermionic action is the Wilson action for the massless 
Schwinger model (k — 1/4). 

1. We generated 50 gauge field configurations A' on the coarse 7x7 lat- 
tice according to their probability distribution e _/3 Sp ( A ' (with periodic 
boundary conditions). For definiteness this has to be done in a fixed 
gauge. Technically, we do this by randomly sampling the non-zero di- 
agonal elements of the momentum space representation of the action 
[21]. For each of these configurations we then calculate the correspond- 
ing minimizing configuration A m i n (A') as described above. 

2. With the help of A m [ n (A') we generate the fermion matrix on the fine 
lattice and performed the BST (Grassmann integral) for the fermions 
giving the (2 • 7 2 ) x (2 • 7 2 ) fermion matrix MBST{A m i n ) on the coarse 
lattice. This is done for all 50 gauge field configurations. 

3. The resulting fermion matrices are then compared with the fermion 
matrices Mp(A') for the coarse lattice generated for the correspond- 
ing A' configurations. A new set of parameters according ()16|) is now 
determined by minimizing 



J2\\ m bst (A min (A')) - M F (A')\\ 2 , 



(19) 



A' 



where the matrix norm is defined 




(20) 



10 



All these steps are iterated until the coupling constants remain stable within 
small statistical fluctuations. 

We worked at f3' = 20 which corresponds to a typical value for the field 
strength F(x) ~ 0.18. After O(10) iterations the coupling constants sta- 
bilized and subsequent iterations stayed within the statistical errors. The 
largest matrix element is ~ 1.4; the average deviation between individual 
matrix elements given by the BST and by the parameterization was 0.00024 
in the last iteration step. Comparing different sets of A' the values of the 
coupling constants we obtained in this way differed in the mean by 0.0017, 
which provides an estimate of the possible systematic uncertainty due to the 
finite number of configurations. We have no control on possible redundancies 
in the parameterized action. Thus it may well be, that part of the observed 
(small) fluctuations in the couplings is due to cancellations of certain terms 
in the fermionic action. The final action therefore may have been determined 
to an even higher accuracy. 





Figure 2: (a) Symbolic representation of some of the 03 couplings of the 
perfect action. Crosses and circles mark the position of the \l/ and Each 
block of arrows corresponds to a term ps(f) 03 U(x, f) ty(x + 5f), where 
/ is symbolized by arrows and the coupling ps(f) is written at the height of 
the arrows, (b) Terms building the perfect clover term. 

In figure El we show a subset of the couplings of the FP-action, which is 
proportional to \j/cr 3 \l/. In the limit of small A^(x) one can expand U(x,f) 
and the expressions linear in A^(x) can be compared with Fig. 14 of ^7]. If 
we compare the two most important operators of this perfect clover term, our 
result agrees within a few percent with • The smaller coupling constants 
of the clover term have already rather large statistical errors. The locality of 
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our FP-action is established by comparing our couplings with the situation 
of vanishing gauge fields. We find that our values are similar or smaller than 
those obtained for the free fermion perfect action jH]. The complete set of 
couplings may be retrieved from |2Uj . 

3 Simulation of the FP-action 

For notational simplicity we call now the obtained parameterization of the 
FP-action just the FP-action. For (3 — > oo the FP-action is on the RT by 
construction. For finite values of /3 it may be no longer perfect. This depends 
on the system considered; the FP-action of the 0(3)-model turned out to stay 
close to perfectness In order to check the amount of improvement, we 
have to rely on direct simulations with the FP-action. Indeed, as will be 
shown below, we found significant improvement for all observables studied. 
One has to determine the path integral 

^ J DA j D^D^ C>(*, tt, A ) e -PS P (A)+S(%%A) _ (21) 

A full-scale (hybrid-) Monte-Carlo simulation for the FP-action is a highly 
non-trivial task. We decided to work on lattices of moderate size up to 
16 x 16. In this case one may perform the Grassmann integrals explicitly, i.e. 
by computing the corresponding determinant and inverse fermion matrix. 

For a simulation we generated 5000 gauge field configurations with peri- 
odic boundary conditions according to their gaussian probability distribution 
e -/3S P (A) with gauge fixing as has been discussed above in Sec. I2.2l in the con- 
text of the BST. For these configurations we determined the fermionic path 
integral as discussed with standard routines of linear algebra. 

The side-benefit of this approach is that we obtain results for both, the 
determinant corresponding to the 1-flavour model and its square, correspond- 
ing to the 2-flavour situation. All observables are obtained in that way. In 
order to estimate the statistical errors we repeated the whole procedure sev- 
eral times (several runs of 5000 gauge configurations each). We compared 
the results obtained with the FP-action with results obtained in in the same 
way for the Wilson-action for fermions (but still using the non-compact form 
of the gauge field action). For the Wilson action we use k = 0.25 through- 
out. This is below the critical value k c (/3) which, however, approaches 0.25 
asymptotically. We therefore do not expect perfect agreement of the obtained 
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physical mass values for thw Wilson action. However, for a qualitative com- 
parison of improvement properties this choice should be sufficient. 
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Figure 3: Improvement of rotational invariancefor the correlation function 
(122(1 in the 1-flavour model for (left-hand plot) Wilson and (right-hand plot) 
the FP-action on a 16 x 16 lattice at (5 = 6. 



Fig. El exhibits for the 1-flavour model the correlation function 

P(x) = (*(0) <r 3 tf(0) V{x) a 3 *(ar)) (22) 

measured for all 2-point separations. The numbers have been determined at 
(3 = 6. One finds much better rotational invariance for the FP-action than 
for the Wilson-action. 

In the 2-flavour Schwinger model one expects one massive mode (which 
we call rj by analogy) and a massless flavour-triplet (called it). The corre- 
sponding momentum-projected operators are 

r](p,t) = y~] e ipXl (u(xi,t) at u(x 1: t) + d{x u t) a x d{x u t)) , (23) 
7r(p,t) = ^ e ipXl (u(x 1 ,t)a 1 u(x 1 ,t) - d(x 1 ,t)a 1 d(x 1 ,t)) . (24) 

Their correlation functions define by their exponential decay the correspond- 
ing energy functions E(p). In Fig0]we present results for the r\ dispersion 
relation for the Wilson-action and the FP-action and compare these with the 
continuum dispersion relation. Again we find significant improvement. 
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Figure 4: Dispersion relation of the massive boson in the 1-flavour Schwinger 
model with Wilson fermions (diamonds) and with the FP-action (full circles) 
calculated on a 16 x 16 lattice at (5 — 6, compared with the continuum 
dispersion relation E(p) = \Jm? + p 2 . 
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Figure 5: Dispersion relation of the massless (tt) boson in the 2-flavour 
Schwinger model with the FP-action calculated on a 16 x 16 lattice at /3 — 6, 
compared with the continuum dispersion relation E{p) = p. 



14 



m 



0.1 



0.5 



0.3 




-0.1 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



17p 



1/2 



Figure 6: Masses for the massive (full circles) and the massless (diamonds) 
modes in the 2-flavour Schwinger model vs. l/vT^i the values have been 
determined with the FP-action on an 16 x 16 lattice. The line denotes the 
theoretical expectations for scaling. 



The dispersion relation for the 7r-state is plotted in Fig. El only for the 
FP-action. For the Wilson action at k — 0.25 and (3 = 6 this state has a 
non- vanishing mass and is therefore not suitable for a comparison. 

It is remarkable that the 7r-mass obtained by the means of the FP-action 
is even numerically very close to zero (e.g. 0.005 at (3 = 6), whereas with 
Wilson-fermions one has to search for a critical point in k. In fact for (3 — > oo 
this constitutes a test for the perfectness of the FP-action, since for the op- 
timal action of course one has to recover the massless continuum Schwinger 
model properties. Also the masses obtained for the massive mode are within 
the statistical errors in agreement with the theoretical expectations, e.g. 
m v = 0.33(1) compared to 0.3257 from continuum theory (at (3 = 1/g 2 = 6). 
Since the mass from the lattice propagator has been determined by the usual 
cosh-fit on 16 x 16 lattices the value has to be handled with some caution. 

Fig. |U1 shows the obtained masses for massive and massless modes at 
various values of (3 (determined with the FP-action on 16 x 16-lattices). We 
find clear signals of non-vanishing 7r-masses at sufficiently small (3, indicating 
deviation of our FP-action from the renormalized trajectory. However, the 
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overall scaling behaviour predicted (for the 2-flavour model) from theory, 



is nicely recovered for moderately large values of (3 > 3. 

4 Conclusion 

For the 2D Schwinger model we determined the optimal FP-action for in- 
teracting fermions in a non-perturbative background of gauge fields. The 
parameterized action contains 123 independent couplings localized in a 7 x 7 
lattice centered around one of the fermionic fields. 

With this action we have simulated the model on a 16 x 16 lattice for 
various values of the gauge coupling and for one and two species of fermions. 
Our results show excellent improvement of the important propagator observ- 
ables, although we did not improve those operators. 

Acknowledgment: We wish to thank W. Bietenholz, F. Farchioni, I. Hip, 
U.-J. Wiese and E. Seiler for discussions. 
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